Free expansion of lowest Landau level states of trapped atoms: 
a wavefunction microscope 
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We show that for any lowest-Landau-level state of a trapped, rotating, interacting Bose gas, the 
particle distribution in coordinate space in a free expansion (time of flight) experiment is related to 
that in the trap at the time it is turned off by a simple rescaling and rotation. When the lowest- 
Landau-level approximation is valid, interactions can be neglected during the expansion, even when 
they play an essential role in the ground state when the trap is present. The correlations in the 
density in a single snapshot can be used to obtain information about the fluid, such as whether a 
transition to a quantum Hall state has occurred. 



There has recently been interest in states of trapped 
atoms that are rotating near the critical frequency of a 
harmonic trap [1—9]. When the atoms are bosons, there 
is a vortex lattice at moderate rotation rates [2], but at 
high rotation rates quantum fluctuations are predicted to 
melt the lattice [5,8] and produce exotic highly-correlated 
many-body states that are related to those occurring for 
electrons in the quantum Hall effect [1,4-7]. A standard 
experimental technique is to turn off the trap potential 
suddenly, and then take a snapshot of the cloud of atoms 
after it has expanded to many times its original size (a 
free-expansion or time-of- flight experiment). A question 
that arises is what information can be extracted from 
such an image of a quantum Hall state. In this note, 
we point out that in a relevant regime, the density dis- 
tribution in such a snapshot directly represents the den- 
sity at the time of switch-off, but rescaled and rotated 
by 90° about the original rotation axis. This result is 
simpler than the results of previous work on the free ex- 
pansion of a Bose condensate in the Gross-Pitaevskii or 
Thomas-Fermi regimes, that can describe a Bose con- 
densate or vortex lattice. While in a few cases the length 
scales of the density distribution are simply scaled up 
by a factor, more generally the condensate function also 
evolves during expansion [10-13]. Usually, the form of 
the density distribution after a free expansion is not sim- 
ply related to the initial density; rather, for noninteract- 
ing particles it is related to the momentum distribution. 
Although our result can be extracted from special cases 
of earlier analyses, the emphasis is different. We em- 
phasize that the neglect of interactions during the free 
expansion is justified whenever all the bosons can be as- 
sumed to be in the lowest Landau level (LLL) just before 
the switch-off, even though interactions may be crucial 
in the highly-correlated ground state that exists before 
that time. Consequently, the result applies to any LLL 
many-body state, not only a Bose condensate. (These 
points were also made briefly in Ref. [7] for the density 
expectation value.) We also suggest how a snapshot of 
the particles after a free expansion can be used to ob- 
tain information about the nature of the original state, 
as it constitutes a microscope that enlarges the real-space 



image of the state. 

We begin by explaining the classical version of the 
problem, which will make the later quantum many-body 
treatment easily understandable. We take particles of 
mass M in a harmonic trap, so the particles move in 
a three-dimensional harmonic oscillator potential, with 
frequency L03 for oscillation in the 3 direction, and u>± 
for oscillation in the 1-2 plane. Thus in general our 
model has rotational symmetry about the 3 axis only. 
Consider a single particle in this potential. The gen- 
eral form of an orbit is an ellipse centered at the center 
of the trap. However, we wish to focus on the single- 
particle states referred to in the quantum case as the 
LLL. These are the states of lowest kinetic energy for 
each positive value of the angular momentum about the 
3 axis. The corresponding classical orbits are circles in 
the 1-2 plane, circling the 3 axis in the positive direction. 
For such an orbit, the momentum when the particle is at 
r = (xx,X2, X3) is clearly p = Mu>±i\j, x r, where A3 is a 
unit vector in the 3 direction, and the angular momentum 
is r x p = Mw^|r| 2 ri3. If the trap potential is removed 
at time t = 0, at which r = ro, p = Po, the particle 
travels freely in a straight line. At time t it will be at 
ro + pot/M = ro + uj±th.3 x ro ~ oj±tn.3 x ro for large 
t. Hence if we have a collection of particles all moving 
on such orbits, their distribution in the 1-2 plane a long 
time t > 1 after removing the trap potential will be 
the same as it was at t — 0, except for a rescaling by u±_t, 
and a rotation by 7r/2 about the 3 axis. In this argument 
we have neglected the interactions between the particles 
once the trap potential is turned off. We will return to 
this question after discussing the quantum case. 

We now turn to the fully quantum-mechanical treat- 
ment, beginning again with a single particle. Suppose 
the state at time t — is ip(r, 0) = ipo(r). For the free 
motion, the solution is simple in momentum space. We 
define the Fourier representation of the state ip by 



i>(r,t) = 
and similarly for ip^. 
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The normalized single-particle basis states for the LLL 
have wavefunctions 



u m (z,x 3 ) 



zme -\z\ 2 /(2l 2 ± )-x 2 3 /(2ll) 
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Here 2: = xi+ix2, the angular momentum quantum num- 
ber to = 0, 1, 2, ... , and l± = y/h/Muj^, l 3 = ^Jh/MLo 3 . 
We note that these states have their maximum amplitude 
on the circles in the 1-2 plane of radius \z\ = l±y/m, 
which is the same as that of the circular classical orbit 
with angular momentum rah. To calculate the Fourier 
transforms, it is useful to introduce the generating func- 
tion 



h {z,xz)=e^ l --\>\ 2 l^-**/^\ 



(4) 



of which the mth derivative with respect to £ is propor- 
tional to u m (z,x 3 ). The Fourier transform of ho is 



(27r) 3/2 /i/ 3 e~^ i - |fc|2 'i /2 - fc 3 i 3/ 2 , 
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where k = k\ + iki- Using (2), at time t, the generating 
function evolves to 



h(z,x 3 ,t) = 

e £z/[l ± (l+i 



\2i]_(l+i 



(1 +iu±t)(l + ioj 3 t) 1 / 2 
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It follows that if a general normalized LLL initial state 
is written as tp (z,x 3 ) = f{z)e~\ z \ /( 2l ±)- x 3/( 2l 3) ; where 
f(z) is complex analytic in z (for example, / can be a 
polynomial), then at time t it becomes 



ip(z,x 3 ,t) = 



'■/[2l 2 , (1+i 
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The probability density is therefore 

\iP(z,x 3 ,t)\ 2 = 
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This is simply the initial probability density \ipo\ 2 , with 
a rescaling and a rotation in the 1-2 coordinates, and x 3 
also rescaled. Notice that the factor 1 + iw±t describes 
the same rescaling and rotation of the initial complex po- 
sition zo — z/(l + ioj±t) in the 1-2 plane as the Cartesian 
formula r = r + w±_t n 3 x r that appeared in the clas- 
sical argument. We should point out that our result is 



a special limiting case of expressions obtained by other 
methods [10,12], on ignoring interactions and restricting 
to the LLL, and is implicit in [7]. 

For large t (t> I/W3), we will be interested in 

r of order t, so define V = r/t. As t — > 00 with V fixed, 
ip(r,t) itself is asymptotic to 



i ^-| z | 2 [(i^_ L t)- 1 +(^_ L t)- 2 ]/(2/i) 

iu>±t(iu) 3 t) 1 / 2 
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In this limit, the probability density (8) becomes 



(LU ± t) 2 (t0 3 t) 



(9) 



(10) 



The transformations have simplified to a rescaling by tu±t 
and a rotation by 7r/2 in the 1-2 plane, and a rescaling 
by u) 3 t in the 3 coordinate. 

The long-time result can also be obtained by an appli- 
cation of the stationary phase approximation to eq. (2), 
which gives the limit t — > 00 with V = r/t fixed, which 
for a general state tpo is 

, (M) ^(f) e ™-^(_iL) s/2 , (11) 



and for the probability density 
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These well-known expressions mean that at long times, 
the distribution in position space is determined by the 
initial momentum distribution, as if the particle propa- 
gated classically with velocity V = hk/M. This result is 
familiar in optics as the Fraunhofer limit of diffraction, 
for example in a two-slit experiment. In the LLL case, 
using eq. (5), the Fourier transform of the LLL initial 
state is 



^o(k) = {2irfl 2 llhf{-ikll)e 
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This shows that for the LLL, the probability density in k 
space is the same as that in position space, rescaled and 
rotated by tt/2. This corresponds both to the classical 
argument above, and to the fact that in the LLL, the 1— 
2 coordinates are canonically conjugate, so x\ = p2l 2 _/h, 
X2 = —pil^/h. Using this in the stationary phase for- 
mula, we recover eq. (10). 

Both the exact and approximate treatments generalize 
trivially to the case of N non-interacting particles with 
an initial state Vo( r ii • • • , r iv)i when it is a linear com- 
bination of products of LLL single-particle states. This 
gives the joint probability density for all N particles af- 
ter the free expansion. This result also holds when the 
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initial state is described by a density matrix within the 
LLL, such as a state of thermal equilibrium, but we will 
not consider that in detail here. 

As a first example, consider a Bose condensate within 
the LLL. Its wavefunction at time t = is a product, 

N 

JlM*], *aj) ( 14 ) 
.7=1 

for some function tpQ. In the LLL, if / in ipo is a poly- 
nomial, one can always factorize ipo{z, £3) oc Yi a ( z ~ 
w a )e _ ' z ' 2 / (2 'i^ a: 3^ 2i ^, and then we have N v vortices at 
complex positions w a , a = 1, . . . , N v . After a long free 
expansion, the wavefunction has the same product form, 
with the bosons condensed in 

l[[z - ;(u,^ a ]e^i 2 [^')~ 1+ ^')~ 2 ]/( 2 <i) 

a 

Thus the vortices will be clearly visible in an image of 
the density in position space, provided the number of 
particles per vortex is large (see the discussion of such 
images below) . This is the condition of large filling factor 
[5] , which is the regime in which a condensed ground state 
of the type (14) occurs [5,8]. 

As a second example, we consider the Laughlin state 
I 14 ], n r<s (*r - zs) 2 ■ n r e-W 2 /( 2 *i)-^/( 2 '3\ where r, 
s = 1, . . . , N. The probability density after a long time 
t is the same as at t = 0, up to the usual two rescalings 
(the rotation can be ignored since the Laughlin state is 
rotationally invariant). This implies that all the spatial 
correlations are preserved after free expansion for time t. 

In our treatment we have completely neglected interac- 
tions during the expansion. As justification for this, we 
note that, when the trap potential is present, the LLL 
approximation should be valid when the s-wave scatter- 
ing length a and typical number density n in the bulk 
of the fluid satisfy 4Trh 2 afi/M < 2hw±, hw^. This condi- 
tion ensures that corrections due to quantum-mechanical 
mixing of non-LLL states into the ground state can be 
neglected. Because n is reduced by centrifugal effects as 
the rotation rate increases, this condition is more likely 
to be satisfied at larger values of the total angular mo- 
menta than at small. (For quantum Hall states, the filling 
factor v — ir 3 '/ 2 l 2 L l^n must be of order around 10 or less 
[5,8].) Even though the interactions are weak compared 
with the kinetic energies in the trap (which are typically 
many times hiu±), the ground states may be highly corre- 
lated. That is because, if interactions are neglected, there 
are many degenerate states of the same total angular mo- 
mentum [1]. (Classically, the particles in the "LLL" or- 
bits keep the same relative positions, up to a rotation, as 
they move, and so the interactions are resonant — over a 
long time they can have a large effect.) The interaction 
energy scale is Aith 2 afi/M ', and when the trap potential is 



removed, n decreases in the same way as the probability 
density, eq. (8). Hence the interaction strength goes to 
zero after a time of order max(l/u;_i_, 1/^3), and the cor- 
rection to the wavefunction, 5tp/?p, due to the interaction 
term a long time after switch-off is small compared with 
1 if 47rS 2 an/Mmax(l/o)^, l/o> 3 ) < h. This condition is 
essentially the same as that for the LLL approximation. 
Further, if the interactions are enhanced by the use of a 
Feshbach resonance (but still weak enough to use the LLL 
approx), then they can be greatly reduced at the same 
time that the trap is turned off, improving the accuracy 
of our neglect of interactions during the expansion. 

Finally, now that we have shown how the free expan- 
sion of a LLL state acts as a wavefunction microscope, 
we point out that even a single high-resolution snapshot 
taken after a free expansion of a highly-correlated state 
contains a lot of information that can be used as a di- 
agnostic tool for the many-body physics. Such a snap- 
shot in principle gives a single typical configuration of all 
the particle positions. In view of the preceding discus- 
sion, we can discuss this in terms of the configuration at 
t = 0, which is drawn from the joint probability density 
\ipo(zi, . . . , zjv)| 2 [we will assume the coordinates are pro- 
jected to the 1-2 plane, so we drop X3 and work in two 
dimensions (2D) from here on]. In a highly correlated, 
incompressible state such as Laughlin's, long-wavelength 
density fluctuations are suppressed, and this can be seen 
even in a single snapshot, if the particle number is large 
compared with 1. For example, such a snapshot is shown 
in Ref. [15]. It differs markedly from a random configu- 
ration. The correlations can be quantified by construct- 
ing the two-particle correlation function of the snapshot. 
This is just a histogram of the values of — rj for all 
pairs of particles i, j (it will be automatically invari- 
ant under r — > — r). This can be compared with the 
two-particle quantum- (and thermal-) average correla- 
tion function g(r) = (ifi (z)^ (z')ip{z')ip{z)) /n 2 (ip(z) is 
the 2D field operator, and r = \z — z'\) which has been 
calculated (usually in an edgeless geometry) for various 
incompressible fluid ground states. We note that in the 
thermodynamic limit, g(r) is normalized so that it ap- 
proaches 1 asr^oo. The wavefunction has an ergodic- 
ity property that ensures that even a correlation function 
constructed from a single sample (snapshot) reproduces 
the quantum/thermal-average g(r), provided the particle 
number is large. Further, the Fourier transform of ng(r) 
is essentially the static (i.e. instantaneous) structure fac- 
tor s(q) [16]. 

There are two basic results that can be extracted from 
the density or its correlations g(r) measured for a sin- 
gle snapshot. First, if one examines a subregion (say, a 
square) of area A of the fluid and determines the par- 
ticle number N s < N in this region, this number will 
fluctuate as the subregion is moved over a given snap- 
shot, and also from one snapshot to another. If the 
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side of the square subregion is larger than the length 
£ we define below, but small enough that the mean N s 
is <C N, then in an equilibrium state the fluctuations 
AN S will be of order (AN S ) 2 = k B TAdn/dfi, and give 
an estimate of the thermodynamic compressibility dnjd\x 
{[i is the chemical potential), if the temperature T is 
known. Alternatively, in the thermodynamic limit one 
has lim q ^ s(q) = {ksT '/n)dn/ 1 d\i [17]. In the zero- 
temperature limit, dfi/dfj, goes to a finite (either zero or 
nonzero) value (except in the case of exactly-zero inter- 
actions in a Bose gas!), so these fluctuations vanish. For 
a translationally-invariant fluid in the LLL in the ther- 
modynamic limit, there is a LLL-projected version s(q) 
of s(q), which is easily obtained from the latter [16]. At 
T = 0, s(q) vanishes faster than q 2 as q — > 0, and for 
incompressible fluids (those in which lim^^o dh/d/i = 0) 
it goes as q 4 [16]. The latter behavior is analytic, imply- 
ing that g(r) — 1 tends to zero rapidly at large r. At a 
finite temperature, there will be a correlation length £, 
which diverges as T —* 0, and which is defined by the 
property that s(q) will cross over at q ~ from the 
T > behavior at q — > to the T = behavior at larger 
q. Note that in an incompressible fluid such as Laugh- 
lin's, £ diverges as £ - e AE/(±k B T) as T ^ , where AE 
is the quasiparticle-hole excitation energy or "gap" . The 
behavior at q > is the second property to look for 
using snapshots. Thus from these correlation properties, 
it is possible in principle to distinguish an incompressible 
quantum fluid from a thermally-melted vortex lattice. 

In practice, there will be both a resolution function 
convoluted with the particle positions, and the question 
of the accuracy with which the density is measured at 
any point (noise). However, our proposal utilizes long- 
distance correlations where the spatial resolution should 
not be a problem, while the average over positions in 
a single snapshot alleviates the accuracy problem, and 
reduces the noise, which is presumably uncorrelated. 

We would like to contrast our remarks with a proposal 
of Sinova et al. to measure the condensate fraction from 
the density profile [9]. They point out that, because of 
analyticity properties in the LLL, the diagonal density 
matrix (i.e. the quantum expectation (n(r)) of the single- 
particle density) determines the off-diagonal density ma- 
trix from which the condensate fraction can be defined. 
Their formula in our notation is 

|A|n(q)| s el«/ 4 , (16) 

where n(q) is the Fourier transform of (n(r)). The nor- 
malization is such that the maximum possible value of 
/? is (3 = 1, and is attained in a product state, as in eq. 
(14). With the help of examples [9], it becomes clear that 
when the number of vortices N v , say in a vortex lattice 
state, is large, the integral is dominated by large q val- 
ues (up to |q| around which tend to be suppressed 
by multiplication by the Fourier transform of the spatial 



resolution function. Further, even in a snapshot with 
perfect resolution, the particle density would be a sum 
of (^-functions at the particle positions, and would differ 
from the average density (n(r)} because of the presence 
of quantum (and more generally thermal) fluctuations. 
These fluctuations are larger relative to the mean, (n(r)), 
at small values of v, which is precisely the regime of great- 
est interest. They also remain large at large q (s(q) — > 1 
[16]), whereas n(q) — ► [9], and (as for any noise in the 
determination of the density at r) they are enhanced in 
the integral by the factor e' q ' l ±/ 4 . These fluctuations 
cannot be removed by averaging over space without de- 
stroying the large q information that is needed. It will 
be necessary to average over many snapshots to obtain 
the quantum average (n(r)). 

While we were completing this paper, a discussion with 
some overlap with the second part of ours appeared [18]. 
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